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It is shown that, despite strong nonlinearity, entanglement 
of formation of two-qubit state can be measured without prior 
state reconstruction. Collective measurements on small num- 
ber of copies are provided that allow to determine quantum 
concurrence via estimation of only four parameters. It is also 
pointed out that another entanglement measure based on so 
called negativity can also be measured in similar way. The 
result is related to general problem: what kind of information 
can be extracted efficiently from unknown quantum state ? 
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Two-qubit state entanglement is well characterised. It 
is the only case where analytical formula for entangle- 
ment of formation -E/ [0 has been provided [?] with help 
of so called concurrence (see ||^ for review). The simple 
necessary and sufficient conditions for the presence of en- 
tanglement in two qubit state is known |^,^ involving so 
called positive partial transpose (PPT) map. However 
the question is how to detect the presence of entangle- 
ment in unknown state possibly efficiently i. e. with min- 
imal number of estimated parameters. As one knowns 
PPT test is represented by nonphysical operation. Con- 
currence and entanglement of formation are highly non- 
linear functions too. Apparently it seems that they re- 
quire full prior state reconstruction in general. The same 
might be expected to hold for any other entanglement 
measures [^,^ (see also 0). Indeed direct detection of 
entanglement of formation function |^,^ has succeeded 
only for pure two-qubit states and relied on very special 
property of the function in that case. Even then, how- 
ever, the procedure requires estimation of more than one 
observable ^. 

Quite recently, using formula of the best structural 
physical approximations of hermitian nonphysical maps 
PH , it has been shown |l2j how to detect violation of 
PPT separability test (or any positive map separability 
test) experimentally without any prior knowledge about 
quantum state (for the analysis of the interferometric 
scheme see Q). 

In this work we provide a protocol detecting concur- 
rence of unknown state by estimation of four parame- 
ters in collective measurements of small (not more than 
eight) number of copies. We also point out how to es- 
timate computable entanglement measure ||^ based on 
negativity in a similar way. For unknown state 

this gives quadratic gain in number of parameters if com- 
pared to quantum tomography. The latter protocol is not 
restricted to qubits but is valid for any d ® dl systems. 



We discuss both presented protocols as far as number 
of involved copies are concerned. In particular to es- 
timate two-qubit entanglement of formation one needs 
slightly more number of copies than in state reconstruc- 
tion scheme. It is however different for "computable" 
measure |^ where not only number of parameters esti- 
mated but also number of pairs involved is less than what 
state reconstruction requires. Subsequently we shall de- 
scribe the protocols and conclude with some discussion. 

Measuring two-qubit entanglement of formation .- 
Let us start with general formula for entanglement 
of formation Ef of given bipartite state g [Q: 
Ef{g) = min J:^p,SiTrA{\^J^){^p^\)) where S{a) - 

—Tr{alog2<y) is von Neumann entropy of state a (counted 
in bits) and the minimum is taken over all ensembles 
{Pi,''Pi} such that J2iPi\'4'i){''Pi\ = Q- For two qubits the 
entanglement of formation has been explicitly calculated 
[01 and it amounts to: 
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where Shannon binary information is h{x) = —xlog2{x) — 
(1 — x)log2{^ — x). The important state function called 
concurrence [?] 
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C{g) = max[\/ Xi — \/ X2 — \/ ^3 — \/ A4, 0] 



involves four real monotonically decreasing numbers {A^} 
which are eigenvalues of (nonhermitian) matrix gg with 



g = S 
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where Pauli matrices act on Alice and Bob qubit respec- 
tively (recall that g is two qubit state) and star stands 
for complex conjugate. Obviously neither C nor Ef is 
measurable in usual quantum mechanical sense i. e. nei- 
ther of them is an observable. Indeed this is forbidden by 
the very fundamental laws of quantum mechanics: both 
functions are nonlinear in state parameters while only 
linear operations can be performed on single copy. Still 
one can try to estimate them in another way - involv- 
ing access to several copies. Such approach succeeded for 
pure two-qubit states based, however, on the very spe- 
cial property of any pure state ie. the dependence of its 
entanglement on eigenvalues of reduced density matrix. 

We would have C and Ef determined for arbitrary 
given g if only we knew the four numbers: Ai, A2, A3, 
A4. They come from rather complicated nonlinear func- 
tion of g. Surprisingly for any unknown g they can be 
physically measured by estimating only four parameters 
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(instead of fifteen ones describing full density matrix) if 
joint measurements are available. The idea of the present 
method is to find (in collective quantum measurements) 
four "moments" that allow us to reconstruct A-s and, 
consequently, concurrence completely. The basic scheme 
of the corresponding protocol is very simple: 

g^g^gi =^ (Mi) = ^A, 

i 

i 

" " " " " ^^g3 =^(M3) = ^(A0^ 
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According to the above, given sample of many copies 
of given g we first divide it into four groups consisting 
of two, four, six and eight copies respectively. Then 
we subject each group to some specific quantum chan- 
nel (i. e. completely positive tracepreserving map) A^, 
k = 1,2,3,4. As a result we get quantum states gk- Fi- 
nally we estimate mean values (Mk) = Tr{gkMk) of four 
observables Mk- They represent the "moments" needed 
to reconstruct values {Ai}. Moreover, as we shall see be- 
low, each (Mk) can be determined via estimation of single 
parameter in binary POVM. 

To see how the above scheme works we need to recall 
the concept of physical approximation of PPT operation 
in lab . Namely the tracepreserving map 
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is physically implementable on any bipartite d® d state 
axY defined on Tlx ® Ti-y. Thus Axy corresponds to 
what in literature is called quantum channel. In formula 
(^ Ix , ly stand for identity matrices corresponding to 
subsystem X, Y respectively. The map I is just iden- 
tity map on subsystem X while T stands for transposi- 
tion on second subsystem Y. The partial transposition 
fx ® Ty simply transposes indices corresponding to the 
second subsystem. The map (^ is called structural phys- 
ical approximation of the nonphysical PPT map 1 ® T 
In particular one can define in similar way such 
approximation for transposition map T. For two qubits 
the latter coincides with previously introduced universal- 
NOT gate (see Jl^). Subsequently detailed description 
of the scheme (|4) will be presented in two steps. 

Step I. Action of channels Ak ■- We define quantum 
channel A^ as being composed of two physical opera- 
tions: (i) the map (||) with d = dk = i'' and spaces Tix 
{Hy) describing all the copies on odd (even) position in 



the k-th group and (ii) the subsequent action of unitary 
operation E on each "even" copy. For instance in the 
second group g® g® g® g {k — 2) the first and the third 
copy corresponds to subsystem X, while the remaining 
two represent subsystem Y. Only copies belonging to Y 
will be subject to subsequent action (ii) of unitary oper- 
ation E. 

Let the channels A^; act on the groups of states accord- 
ing to the scheme (^). As a result we get the following 
four states 
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Up to the maximal noise part, these states can be consid- 
ered as "proportional" to the matrices {g® g)^''. In fact 
they differ from the matrices by shrinking factor at 
the "Bloch" . This property is similar to what happens 
in the theory of approximate cloning or in the case of 
universal NOT gate |l|l. 

Step II. Measurement of moments (Mk) ■- We shall 
apply the approach ||ll|,|l^ providing spectrum state es- 
timation with help of observables constructed from shift 
operators Vn Il4|. They have the following property 



Tr{Vi^n)Ai 
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Combining the above property with the fact that all mo- 
ments 'Yliii^iY ^I'S 0116 can see that the latter are 
just by mean values of the following observables 



Mk 
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when each of them is calculated on state gk (see the 
scheme (^)). Finally, the numbers A^ can be calculated 
uniquely from the moments. This can be shown in anal- 
ogy to the analysis of Ref. [0 where state spectrum es- 
timation was provided. 

Now an important point is that each mean {Mk) can 
be detected (up to rescaling) as a single parameter in 
special binary POVM: following Ref. ||l^ each observable 
Mk can be encoded into some ancilla and then its mean 
value can be reproduced form measurement of elementary 
binary observable (Pauli matrix tj^) jl^. This concludes 
the description of the protocol. 

The presented method is parametrically efficient - it 
requires much less parameters to be measured than state 
reconstruction does (four instead of fifteen). It is not so 
as long as number of copies is concerned: one "round" of 
quantum tomography requires 15 copies while here (see 
(^)) we require 2-1-4-1-6-1-8=20 copies in each round of the 
protocol. However the factor r = • = "number of 
parameters x number of copies" is here superior (r = 
80) than for the state reconstruction schemes (r = 165). 
Even more striking is the fact that while the number 
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of copies consumed in one round of the experiment is 
increased only by 33% (from 15 to 20), the number of 
parameters is decreased almost four times (from 15 to 
4) . It suggests that the optimal method of entanglement 
estimation should exists that consumes the same number 
of copies as state reconstruction but requires less number 
of parameters (5 or 6). Note that the present method 
is useful when we are allowed to use small number of 
apparata with fixed architecture. One can also anticipate 
the existence of trade-off between Vp and Tc involved in 
estimation of any entanglement. The optimisation of this 
trade-off is an interesting open problem. 

Note that so far we have asked "How much entangle- 
ment is in the system?" . The above measurement proto- 
col not only answers but, obviously, solves also the less 
detailed problem: "Is there any entanglement in the sys- 
tem at all ?" . One can however discern the two questions 
above and consider the scenario when the observer knows 
that two qubit entanglement is present in the system but 
has no idea "how much" . Below we shall provide protocol 
for such a scenario basing on the fact |l^ that if 2(g)2 state 
g is entangled then the smallest eigenvalue of the matrix 
7 = T^g^^I^g^^ is proportional to . The modi- 

fied protocol must involve the map (||) (with transposed 
subsystems suitably permuted) to produce the new states 
Yk with Bloch vector "proportional" to (Egi^'^E® g^^)®'' 
(that play the role of states gk form Step 7). Then one 
should estimate eigenvalues of "fk in the same way as in 
Step II and infer C{g) following result of Ref. ||l^. The 
above protocol answers quantitative "how much" ques- 
tion provided that observers has qualitative knowledge 
that entanglement is present. The latter must be estab- 
lished before and to this end physical application ||l^ of 
PPT test can be used. Due to its binary character ("Yes- 
No" answer) it requires less measurement precision than 
the qualitative stage described above. We can use it to 
minimise measurement effort: if PPT test does not reveal 
any entanglement we just abandon the second stage. 

Multilevel systems and computable entanglement mea- 
sure .- For general d <S!) d' systems (they can be called 
multilevel systems) there is no analytical formula for en- 
tanglement of formation. However at least one nontriv- 
ial entanglement measure can be calculated analytically. 
This is so called "computable" entanglement measure [|j : 



Ec{g)^log2\\g^ 



(9) 



where g^^ = [I (g) T]{g) and || • || stands for trace norm 
which for hermitian operator means sum of moduli of its 
eigenvalues. In particular the quantity (| If?"^^ 1 1 ~ 1) intro- 
duced to quantify entanglement ||l^] and after normalisa- 
tion N — (I 1 1 — l)/2 called negativity has been shown 
to be entanglement monotone ]20[ |. The measure (^) 
is weaker than entanglement of formation, it "detects" 
however all free ie. distillable entanglement of bipartite 
systems. To estimate the measure we must experimen- 



tally determine the quantity 
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where {A^} are eigenvalues of partially transposed g sym- 
bolised by g^^ . The above formula is crucial for subse- 
quent analysis. 

We now apply scheme of Ref. (l2| making however full 
use of its output rather than focusing on minimal eigen- 
value. For any bipartite d®d state g the scheme consists 
of application the map (H) which leads to the new quan- 
tum state gi = Kxy{q) with eigenvalues {Aj}. Each sin- 
gle A; one eigenvalue A ■ of g^^ . Indeed they are related 
by the following afhne function: 



A,: = 
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Now the key point is that {A^} can be measured exper- 
imentally via measurement oi rp — d^ — 1 observables 
(see pl] , p^p^ ) ■ It is much less than — 1 required for 
quantum tomography. From results of the measurements 
one can calculate all parameters A^ = {d^ -f l)Ai — d and 
substitute them to the formula 1^ which provides the 
needed value of the corresponding measure. 

The whole scenario can be immediately generalised to 
any d ® d' system with d ^ d! . Only the parameters in 
(^ will change slightly and this can be easily calculated 
with help of structural approximation formula [ll| ]. 

It is interesting that in the present scenario one has 
also gain in number of copies involved Indeed we have 

= 2-h3-H...-f'rf2 ^ (^4 _^ ^2 _2)/2 instead ofr^ = d^'-l. 

Discussion .- In conclusion we have provided the first 
nontrivial protocols of entanglement estimation that re- 
quire no prior knowledge about the states. The first 
presented protocol allows to detect two-qubit entangle- 
ment of formation with help of joint measurements on 
small number of systems. In one round the method re- 
quires 33% of copies more than state reconstruction (20 
instead of 15). However it needs almost four times less 
parameters (4 instead of 15) than quantum tomography 
does. We have also shown how to estimate experimen- 
tally computable measure of entanglement Ec based on 
partial transpose operation. Applying the protocol of 
we pointed out that the measure can be determined 
with help of estimation of only d^ — 1 parameters in- 
stead of d'' — 1 ones required for state estimation. The 
scheme provides also gain in the number of copies con- 
sumed: Tc = (c?^ + d^ — 2)/2 instead of Tc = d'* — 1 due 
to reconstruction schemes. 

It would be interesting to find similar protocol deter- 
mining some parameters of the best separable approxima- 
tion of two-qubit quantum states (for instance entangled 
part of Lewenstein-Sanpera decomposition ||2l|] ). It in- 
volves the question how to find experimentally eigenvec- 
tor corresponding to the least negative eigenvalue of par- 
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tially transposed density matrix q^^ . To get the eigen- 
vector with minimal measurement cost would be partic- 
ularly interesting because it plays crucial role in the only 
universal protocol of two-qubit entanglement distillation 
p2[ . In the above context basic questions naturally arise: 
(i) whether it is possible to estimate eigenvector corre- 
sponding to extreme eigenvalue of given density matrix 
without prior reconstruction of the latter ? (ii) if so, 
how efficiently can it be done ? Any possible solution 
of these problems must take into account fundamental 
restrictions on quantum "nonlinear operations" [ pl]j2^ ]. 

In any case the present results lead naturally to an- 
other interesting question: is it possible to estimate ef- 
ficiently other entanglement measures (like relative en- 
tropy of entanglement [^ ) or - more generally - cor- 
relation measures like the one from Ref. |2j]) based of 
von Neumann entropy ? For pure states it can be done 
as all measures are equivalent to the entropy of the re- 
duced density matrices which at the same time is pro- 
portional to index of correlation I = Sa + Sb — Sab- 
For mixed states the latter serves rather as a measure of 
global (quantum -f classical) correlations [Q. Following 
the present analysis the latter can be efficiently measured 
without state reconstruction in general. At present one 
can not develop the method for the relative entropy of en- 
tanglement or distillable entanglement because existence 
of analytical (or partially analytical) formulas for those 
measures is still an open issue. 

Furthermore from the point of view of general quantum 
information theory very important problem is: what kind 
of information ( whatever it means ) can be extracted from 
unknown quantum state at small measurement cost?. 
Here we have considered extraction of information repre- 
sented by special parameter - entanglement measure. It 
is closely related to the following problem what kind of 
measurement information can be extracted from the state 
in the protocols destroying as little quantum information 
as possible ? Indeed measurement of small number of 
parameters can be always interpreted as an incomplete 
von Neumann measurement which destroys quantum in- 
formation less than the complete measurement does. An 
illustrative example is a detection of spectrum without 
state reconstruction via Young diagrams projections pro- 
posed in Ref. (for alternative spectrum detection see 
|p3| ). The above issues are crucial for some quite fun- 
damental processes like universal quantum compression. 
For example it has been shown ||2^ that quantum source 
can be compressed efficiently if only one knowns sin- 
gle parameter represented by von Neumann entropy of 
the source. The entropy can be estimated without prior 
knowledge about the state: the measurement effort (see 
p3[ , cf . ) is the same as in the estimation of the com- 
putable measure considered above. Finally, we have con- 
sidered one of the basic problems in general entanglement 
theory: "how one can detect the amount of entanglement 
experimentally?" . This question concerns not only en- 



tanglement itself: any answer to it says something about 
properties of information contained in unknown quantum 
states. 
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